Competition between Hidden Spin and Charge Orderings in Stripe Phase 
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The correlation between charge and spin orderings in hole-doped antiferromagnets is studied 
within an effective model of quantum strings fluctuating in an antiferromagnetic background. In 
particular, we perform the direct estimation of the charge and spin long-range-order parameters 
by means of the quantum Monte Carlo simulation. A hidden spin long-range order is found to be 
governed by a competition between the two trends caused by increasing hole mobility: the enhance- 
ment of the two-dimensional spin-spin correlation mediated by hole motions and the reformation of 
a strong stripe order. 
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Since the observation of higli-tcmperature supercon- 
ductivity in a hole-doped antiferromagnetic (AF) insula- 
tor, spin dynamics in hole-doped low-dimensional Heisen- 
berg antiferromagnets has turned into a major subject of 
theoretical and experimental studies. The discovery of an 
anomalous stripe phase within Cu02 planes of the Cu- 
oxide superconductor, La2_2;-i/Ndj,Sra;Cu04, has stim- 
ulated considerable interest in the studies of the corre- 
lation between the charge and spin orderings in hole- 
doped antiferromagnets The experimental findings 
suggest that, in the stripe phase, the dopant-holes are 
segregated into one-dimensional (ID) domain walls which 
separate AF anti-phase spin domains 0, Q . The stripe 
phases have been experimentally observed in a wide 
range of doped-hole concentrations for the lanthanum 
cuprate 0,01 and proposed in the two- dimensior ial ( 2 D ) 



Hubbard-like lH 1^121,0 113 and J models [Uli^. 
It is generally accepted that the static stripe order, e.g., 
observed for La2-a;Sr2,Cu04 with x ~ 1/8, is pinned by 
the lattice modulation in the low-temperature tetragonal 
phase. 

The stripe order with AF anti-phase spin domains has 
also been observed in Ni-oxide La2-a;Sra;Ni04 0, Hsf . 
though the material remains insulating except at very 
high Sr concentrations. There is by now a wide consen- 
sus that the formation of stripe is a generic property of 
hole-doped antiferromagnets [ij. The mean-field theory 
for the Hubbard-Peierls model proposed by Zaanen and 
Littlwood suggests a picture in which holes bind to 
domain walls embedded in an AF background. 

Several years ago Zaanen and coworkers had intro- 
duced an appealing model for stripes acting as anti- 
phase boundaries intervening AF spin domains, which 
captures well the essential features of the interplay be- 
tween spin and charge degrees of freedom. This effec- 
tive model, deduced on basis of the Ogata-Shiba prin- 
ciple is described in terms of fluctuating quantum 
strings. Zaanen's model can be viewed as a variant of 
the bond-alternated spin model with bond randomness 
as mentioned below and is interesting in its own rights 
as a new random quantum spin system. Our main pur- 
pose in the present work is to investigate the correlation 



between the spin and charge orderings and the tendency 
of the spontaneous ordering of stripes within the effec- 
tive model. In particular, we carry out a high precision 
numerical investigation into the charge and spin order- 
ings by direct calculation of the order parameters at zero 
temperature in the thermodynamic limit by means of the 
quantum Monte Carlo (QMC) simulation. We find that 
a hidden spin long-range order (LRO) is governed by a 
competition between the two trends caused by increas- 
ing hole mobility: the enhancement of the 2D spin-spin 
correlation and the reformation of a strong stripe order. 

The motive of the aforementioned model is to describe 
the quantum mechanical motion of hard core particles 
forming extended strings, which fluctuate in an AF back- 
ground. The strings have an overall orientation which is 
aligned with the y-direction. Each quantum string con- 
sists of Ly hard-core particles referred to hereafter as 
'holes'. Each hole moves in the a;-direction and an ef- 
fective exchange interaction exists between the two spins 
neighboring a hole. By regarding the hole as residing 
on the link between two spins neighboring the hole, we 
can reduce our system to a spin-only model defined on a 
squeezed lattice as depicted in Fig. The model thus 
obtained is described by the Hamiltonian il7| 
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+ J{1 - {I - a)n(^^ y))S(^^ yy S^x-i-i,y) (1) 

+ J(l — ^ '>T-(x-l,y)^(x ,y+l)) 

where the summation Y]^ „ runs over all the lattice sites 
on the squeezed lattice and S(^x,y) is the S = 1/2 spin 
operator at site {x,y). The first term corresponds to the 
kinetic energy of holes, where a|^^ is the creation op- 
erator of the hole sitting on the link between the sites 
(x, y) and {x -\- l,y) in the squeezed lattice, and P is the 
projection operator ensuring that strings are not broken 
up and are separated by at least one spin site in the orig- 
inal lattice. The second term describes the 5=1/2 ex- 
change interactions in the x-direction. We fix J to unity 
hereafter. The value of a corresponds to the effective in- 
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FIG. 1: Transformation from the original lattice to the 
squeezed lattice. The quantum strings in the original lattice 
are described by the dashed lines connecting links denoted by 
bold lines in the squeezed lattice. 

teraction between spins neighboring a hole and n^^ y) is 
the hole number operator at site {x,y). The third term 
makes the exchange interaction zero in a link parallel to y 
in the squeezed lattice for bonds connecting to the holes 
in the original lattice. 

Iir the perfect stripe system with t — 0, in which case 
the quantum strings are perfectly straight and line up at 
regular intervals, the system as described on the squeezed 
lattice reduces to a spin ladder model with an inter- 
ladder interaction of strength a When L^/Ng is 
an odd integer, where is the size of the lattice in the 
X direction and Ng is the number of strings, this sys- 
tem exhibits an even-leg-ladder-like behavior: there is a 
quantum phase transition between the disordered phase 
with a finite spin gap for small a and the AF LRO phase 
for large a |20j. The transition point is located at e.g. 
ac = 0.3138(1), 0.0787(2), and 0.0153(1) for the con- 
centrations X = Ns/Lx = 1/3, 1/5, and 1/7, respectively. 
These results are obtained by the finite-size scaling (FSS) 
analysis of our QMC data of the correlation length with 
the exponent i' — 0.71 fixed to the value of the three- 
dimensional classical Heisenberg universality class [2]| . 
This disordered phase persists in the presence of the hop- 
ping interaction up to a certain value of tc as mentioned 
below. On the other hand, when the value of L^/Ns is 
even, the system behaves like an odd-leg spin ladder: the 
system with a = is in a gapless phase and AF LRO is 
induced by an infinitesimal value of a. 

The strings in the present model, subject to quan- 
tum fluctuation viz. the t-term, are at the same time 
in contact with the AF background, and thus induces an 
imaginary-time dependence of the spin exchange coupling 
on the squeezed lattice. To deal with this aspect, the spin 
configuration is therefore updated by the discontinuous 
imaginary-time loop algorithm 22] with fixed hole config- 
uration. Meanwhile, the hole configuration is described 
by the Suzuki- Trotter decomposition via the prescription 
due to Eskes et al. and updated by the Metropolis 
algorithm with a fixed spin configuration. The QMC sim- 
ulation with the loop algorithm combined with the con- 
ventional world-line algorithm is carried out on Lx x Ly 
{Lx = Ly = L < 36) square lattices with the periodic 
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FIG. 2: Plots of (a) the charge-charge structure factor and 
(b) the spin structure factor in the original lattice for x = 1/3 
at T = 0.01 in the charge-ordered and hidden AF long-range 
ordered phase with a = 0.1 and t = 1. Note that the peak at 
(0,0) of the charge-charge structure factor is subtracted. 
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FIG. 3: Plots of (a) the charge-charge structure factor and 
(b) the spin structure factor in the original lattice for x = 1/4 
at T = 0.01 in the a — 0.1 and t = 1 system. 

boundary condition. The value of P/Nt is fixed to be 
1/4, where /3 is the inverse of temperature T and Nr is 
the Trotter number. The error bars are estimated based 
on about 10^ samplings. For each sample, 10^ Monte 
Carlo steps (MCS) are spent for measurement after 10''' 
MCS for thermalization. 

The ground-state phase diagram parameterized by a 
and t for x — 1/3 was previously presented by Zaanen 
et al. 17] . Our result is qualitatively consistent with 
theirs: the system is in the charge-ordered and spin- 
disordered phase for small t and small a and the phase 
transition to the charge-ordered and hidden AF LRO 
phase occurs as the value of t or a is increased. Here- 
after, the AF LRO in the squeezed lattice will also be 
referred to as the hidden AF LRO in terms of the orig- 
inal lattice. The spin-disordered state is effectively the 
same as the state realized in the spin-1/2 even-leg-ladder 
system with the inter-ladder interaction a as mentioned 
above. When a = 0.1 and x = 1/3, our FSS analy- 
sis of the Binder parameter 01 shows that the phase 
transition occurs at tc = 0.21(1). The charge-ordered 
and hidden AF LRO phase for i > tc is characterized by 
peaks of the charge-charge structure factor Sc{kx,ky) = 
X^dSr '^''"''^r]^) at {kx,ky)={±^TT,0), whcrc is the 
hole number operator at site r, and those of the spin 
structure factor Ss{kx,ky) = ^(L-jVs) ^Er e''^''"-5'r]^) at 
{kx,ky)={±^TT,±TT) as shown in Figs. [21(a) and (b). 
The results show that the hidden AF correlations with 
anti-phase boundaries are present, consistent with the 
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FIG. 4: The t-dependence of the stripe LRO parameter at 
T = for a; = i. The square and circle denote the QMC 
results of the stripe LRO parameter for a — 0.1 and 0.5, 
respectively. 



stripe structure, namely, the peaks of Sc{kx,ky) and 
Ss{kx, ky) are at {kx, ky)={zL2Trx, 0) and (±(1 — x)7r, ±7r), 
respectively. Note that there are no additional peaks 
corresponding to other modes, except for the one at 
{kx, ky) — (0,0) for Sc{kx, ky) in the thermodynamic limit. 
While the x = ^ system exhibits a even-leg-ladder-like 
behavior, the system at a; = | is odd-leg-ladder-like. 
This is evident from the peaks that are observed at the 
positions (kx, ky) = {zL^TT,0) and (±|7r, ±7r) as shown in 
Figs. ISl(a) and (b), respectively. It shows that at least 
in the parameter range we calculated, the stripe con- 
figuration consisting of only three-leg-ladders is realized 
rather than an alternating series of two- and four-leg- 
ladders, which was suggested in Ref. where interac- 
tions from holes were neglected. The peaks of Sc{kx, ky) 
at fcy = shows that each string selects a overall direc- 
tion in space, though it can locally fluctuate. This is 
associated with the "directedness" observed in a quan- 
tum string model [i^. On the other hand, the peaks at 
kx = i:2TTx corresponding to strings lined up at regular 
intervals might have their primary origin in the effective 
repulsion coming from the constraint that the stripes be 
separated by at least one spin site. 

We make direct estimations of the f-dependent stripe 
LRO parameter defined by 



(^stripe) 



lim lim 



Sc{2ttx,0) 



(2) 



The value of iS'c(27ra;, 0) converges to its zero-temperature 
value at T lower than an energy scale related to the finite- 
ness of the system size. Thus, 5'c(27ra::, 0) at low tem- 
peratures where its T-dependence becomes indiscernible 
within the error bars is taken as the values in the T ^ 
limit. Furthermore, the value in the limit of L — > oo is 
obtained by fitting the QMC data to ^/ Sci2nx,0)/ L'^ ~ 
(Ostripc) + a/L. The results of the extrapolation are 
shown in Fig. ^ When the system is in the perfect 
stripe phase with t — 0, the value of (Ostripo) equals x, 




FIG. 5: The t-dependence of the hidden AF LRO parameter 
at T = for X = i. The open square and circle denote the 
QMC results of the hidden AF LRO parameter for a = 0.1 
and 0.5, respectively. The filled square denotes the critical 
point tc — 0.21 for a — 0.1. The dashed line is obtained 
by least-squares fitting with the fitting function Mh = c{x — 
0.21)'' in the range of 0.3 < t < 1. The values of c and 6 are 
estimated to be 0.30(1) and 0.23(4), respectively. 

because the strings are perfectly straight and align at reg- 
ular intervals, from which follows that S'c(27ra::, 0) = N^. 
For t < 3 the value of (Ostripc) reduces as t is increased. 
This is consistent with the picture that hole fluctuation 
competes with the stripe order. For t > 3, however, 
the strong stripe order is reformed and the ground state 
approaches the perfect stripe. This reformation of the 
strong stripe will be discussed later. 

We also make direct estimations of the ^-dependent 
hidden AF LRO parameter defined by 



lim lim - 



/ 3S'(7r,7r) 



(3) 



where S{Tr, tt) is the size-dependent staggered structure 
factor in the squeezed lattice. The extrapolation is per- 
formed by a procedure similar to that for (Ostripo)- A 
finite Mh is induced at t > tc = 0.21(1) for a — 0.1, with 
the value increasing with t in the region tc < t < 3 a,s 
shown in Fig. [S] This result suggests that the hole mo- 
tion enhances the effective spin-spin interactions, which 
in turn obscures the ladder-like structure, and yields in- 
stead the 2D magnetic LRO. To paraphrase, the mag- 
netization which had been suppressed in the underlying 
bond-alternating spin system due to quantum fluctua- 
tions is recovered by the randomness of the spin interac- 
tion brought on by turning on a finite t. Randomness- 
induced-LRO is expected to occur also in 2D quantum 
spin systems with bond randomness |25| . The present 
case differs in that the randomness is present also in 
the imaginary-time direction, i.e., one may view it as 
a "2-l-lD random quantum spin system". A further in- 
crease in the hole motions, however, eventually leads to 
the reduction of the spin LRO due to the reformation of 
the strong stripe order, leading to the spin system with 
alternating structure which was seen for < t < tc- The 
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FIG. 6: The t-dependence of the LRO parameters at T = 
for a; = i. The square and circle denote the values of the 
stripe LRO and the hidden AF LRO parameters for a — 0.1. 



properties mentioned above are also seen for a = 0.5, 
where the system with t — > does not belong to the 
spin-disordered phase. 

The values of (Ostripo) for a = 0.1 are larger than those 
for a = 0.5 as shown in Fig. 01 This result is in conflict 
with the following naive energetics: a kink, arising when 
one hole in the straight string moves to a neighboring site, 
yields an energy loss proportional to a nearest spin-spin 
correlation. In the hidden AF LRO phase, a naive mcan- 
field-like argument gives the correlation proportional to 
My. This implies then that in order to gain energy the 
value of the stripe LRO parameter for a = 0.5 should 
become larger than that for a = 0.1 because of the a- 
dependence of Afn. Our results indicate that the specula- 
tion is too naive to explain quantitatively the magnitude 
of the order parameters. 

The i-dependences of (Ostripc) and Afn for a = 0.1 
and X — J are also shown in Fig. El where the sys- 
tem in the limit of i ^ becomes odd-leg-ladder-like. 
Though a similar upturn behavior of the order param- 
eters to that for a; = 1/3 is observed, the i-dependence 
is indiscernible due to the weak effective repulsion be- 
tween quantum strings. We thus expect in general, that 
it becomes harder in practice to see the correlation be- 
tween spin and charge orderings as one goes to lower hole 
densities, in particular for the odd-leg-ladder like system. 

To conclude, we have demonstrated that the magni- 
tude of the hidden AF LRO is governed by the competi- 
tion between the enhancement of the 2D spin-spin corre- 
lation and the reformation of the strong stripe order. The 
question as to why the stripe tends to reorder at large t 
remains to be resolved. Although we have not yet found 
a rigorous description of the reformation of the strong 
stripe order, it may be due to the effective Coulomb re- 
pulsion coming from the constraint that the stripes are 
separated by at least one spin site. The present model 
with the constraint is not realistic at large t. There are 
two further issues along the line of the present study left 
for future work. One is the effect of adding a next-nearest 



neighbor interaction on the squeezed lattice. In that case, 
the diagonal stripe, observed in Cu-oxides and Ni-oxides, 
would be expected to emerge in a certain parameter re- 
gion. The other is to incorporate the dislocation of the 
strings. In the present work, we have restricted the string 
mobility so as to ensure that the strings are not broken 
up. In order to access a superconducting regime, how- 
ever, it would be crucial to study the situation where the 
stripes start to get destroyed [l7]. The numerical simula- 
tion for this situation is a major challenge and new ideas 
are in need. 
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